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Abstract 


In this note, Timoshenko beams with interior damping and boundary damping are studied from the 
viewpoints of control theory and numerical approximation. Especially, the uniform exponential stabil- 
ities of the beams are studied. The meaning of uniform exponential stability in this paper is two-fold: 
The first one is in the classical sense and also is concisely called exponential stability by many au- 
thors; The second one is that the semi-discretization systems, which are derived from an exponentially 
stable continuous beam by some semi-discretization schemes, are uniformly exponentially stable with 
respect to the discretized parameter. To investigate uniform exponential stability of continuous and 
discrete systems, five completely different methods, which are stability theory of port-Hamiltonian 
system, direct method of Lyapunov functional, perturbation theory of Co-semigroup, spectral analysis 
of unbounded operator and frequency standard of exponential stability for contractive semigroup, are 
involved. Especially, a new method, which is based on the frequency domain characteristics of uniform 
exponential stability of Co-semigroup of contractions, is established to verify the uniform exponential 
stability of semi-discretization systems derived from coupled system. The effectiveness of the numerical 
approximating algorithms is verified by numerical simulations. 


Keywords: Timoshenko beam, exponential stability, semi-discretization, finite difference, 
Co-semigroup. 


1. Introduction 


In this paper, we study classical Timoshenko beams with interior damping and boundary damping, 
which was extensively investigated in the past decades due to its widely applications in engineering 
[1-12], both from PDEs and numerical approximation point of view. Adding the effect of shear as well 
as the effect of rotation to the Euler-Bernoulli beam, Timoshenko developed a beam model which is 
depicted by a family of equations as follows: 


{ pwyu(s,t) — K(w,(s,t) — o(s,t)), =0, s € (0,1), t>0, 


Teduti,t)= los, (4,0)— eu, (62) — de.4) = 0, (1.1) 


where t is the time variable and s is the space coordinate along the beam in its equilibrium position. 
w(s,t) is the transverse displacement of the beam and ¢(s,t) is the rotation angle of a filament of the 
beam. The coefficients p, Ip, Æ, I and K are the mass per unit length, the rotary moment of inertia 
of a cross section, the Youngs modulus of elasticity, the moment of inertia of a cross section and the 
shear modulus, respectively. 
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Many authors tried to stabilize the Timoshenko beam by designing boundary damping or interior 
damping in the past decades. Kim and Renardy [1] imposed the boundary conditions 


w(0,t) =0, 6(0,t) =0 (1.2) 
at s = 0 and the boundary damping 
wa(1, 6) H o(1,t) = —aıpw(1, t), ¢s(1, t) E —az2lpġ:(1, t). (1.3) 


at s = 1 for aj, a2 > 0. These boundary damping corresponds to a feedback mechanism which 
monitors w+(s,t) and ¢;(s,¢) at s = 1 and transforms them into the lateral force and moment applied 
at s = 1, respectively. They showed that the energy decays exponentially by the semigroup theory 
and the method of Lyapunov functional. Xu and Feng [2] gave detailed spectral analysis and Riesz 
basis property of the generalized eigenvector system of the beam (1.1)-(1.3). Villegas, Zwart, Gorrec, 
and Jacob et al. [3, 6] obtained the exponential stability of (1.1)-(1.3) by the stability theory of port- 
Hamiltonian system on infinite-dimensional space. Raposo et al. [7] added frictional dissipative terms 
wz(s,t) and (s, t) into (1.1) to get 


pwu(s,t) — K(wes(s,t) — (s,t))s + wi(s,t) =0, s € (0,1), t> 0, 
I, du(s,t) — El dbss(s,t) — K(ws(s,t) — o(s,t)) + (s, t) = 0, (1.4) 
w(0, t) = ws(1,t) = (1, t) = s(1,t) = (0, t) =0, 


which is a Timoshenko beam with interior damping, and obtained the exponential stability of (1.4) by 
the frequency domain characteristics of exponential stability of Co-semigroup. Yan et al. [8] studied the 
stabilization problem of Timoshenko beam (1.1) in the presence of linear dissipative boundary feedback 
controls and presented various necessary and sufficient conditions for the system to be asymptotically 
stable. The equivalence between the exponential stability and the asymptotic stability for the closed- 
loop system was finally given by the same method of [7]. 

Recently, Rivera and Naso [9] imposed boundary dissipation only on one side of the bending 
moment, i.e., HI¢,(0,t) = —a2¢:(0,t), and other three Dirichlet boundary conditions to (1.1). They 
proved the exponential stability of the resulting Timoshenko beam, provided the wave speeds of the 
system are equal. Almeida Júnior, Ramos and Freitas [10] found some new facts related to the classical 
Timoshenko system. More precisely, they proved the damped shear beam model, which corresponds to 
a part of the classical Timoshenko beam model, possessed an energy exponential decay under boundary 
conditions of Dirichlet-Neumann. 

In this note we shall study the beam (1.1)-(1.3) and (1.4) in both continuous and discrete levels. 
In the continuous case, the exponential stability of (1.1)-(1.3) is analyzed by the stability theory of 
infinite-dimensional port-Hamiltonian systems [3, 6]. However, the exponential stability of (1.4) is 
investigated combining the spectral analysis (see e.g. [11, 12]) and perturbation theory of exponential 
stability Co-semigroup [13]. It should be pointed that the results of exponential stabilities of (1.1)-(1.3) 
and (1.4) are not new, we include them to compare with the discrete results. 

Whereas in the discrete case, in which main contributions of this paper are presented, a spacial 
semi-discretization process is firstly carried out for (1.1)-(1.3) or (1.4). A family of discrete systems 
are obtained and a natural problem is posed: whether or not they have an exponential decay rate 
independent of the step size. This is an important problem in studying of dynamical system, see 
for instance [14, 15] and the references therein. To answer this question, we should overcome two 
difficulties. The first one is how to choose an appropriate algorithm from many discretization methods. 
The second one is how to verify that the discrete systems possess uniform exponential stability for a 
given discrete scheme. Fortunately, some related researches [16-22] have been carried out for wave 
equation since 1990s and a series of successful results have been derived. For example, a new kind of 
central difference [20-22] was introduced to bypass the first obstacle and the multiplies method was 
widely and conveniently used to verify that the discrete systems possess uniform exponential decay. 


However, a relatively small amount of work has been made for the Timoshenko beam system, see [23- 
25] for uniform controllability, [26] for uniform observability, and [27] for uniform exponential stability 
of Euler-Bernoulli beam. 

We shall fill the gap and investigate the uniform exponential stability of Timoshenko beams (1.1)- 
(1.3) and (1.4). For this purpose, the Timoshenko beam (1.1)-(1.3) and (1.4) are reduced to a family 
of first order hyperbolic PDEs and they can be rewritten as the form of port-Hamiltonian system 
on a Hilbert space [6, Section 7.1]. A new discretized scheme, which is called the average central- 
difference, is introduced to the order reduction systems motivated by [20-22, 28]. However, to the best 
of our knowledge, suitable discrete multipliers are hard to find for the discrete systems because the 
Timoshenko beam (1.1) is not a simple combination of two wave equations but contains coupled terms. 
Therefore, applying the results of the wave in [20-22] to Timoshenko beam is nontrivial and we should 
deal with it with additional efforts. Finally, perturbation theory of semigroup and multiplier method 
are combined to overcome the second challenge for the semi-discretization systems of (1.1)-(1.3). But 
a new method, which is based on the frequency domain characteristics of uniform exponential stability 
of Co-semigroup of contractions (see e.g. [29, P.162] and [30]), is established to verify the uniform 
exponential stability of semi-discretization systems of (1.4). 

The structure of this paper is organized as follows. In the Section 2, the Timoshenko beam (1.1)- 
(1.3) is studied in both continuous and discrete case. In the Section 3, the Timoshenko beam (1.4) is 
studied similarly by different methods. In section 4, we give some concluding remarks. 


2. Related results of the beam (1.1)-(1.3) 


We discuss uniform exponential stability of Timoshenko beam (1.1) associated with the boundary 
conditions (1.2)-(1.3) in this section. Main results related to the continuous case and its discrete case 
are given in Subsection 2.1 and Subsection 2.2, respectively. 


2.1. Exponential stability of continuous beam (1.1)-(1.8) 


To investigate uniform exponential stability conveniently in the next subsection, we take the trans- 
formation 


yi(s, t) = ws(s, t) = o(s, t), y2(s,t) = pwz(s,t), y3(s, t) = os(s, t), ya(s, t) a pls, t), 


and obtain equivalent system of (1.1)-(1.3) 


jı (s, t) = B2y5(s, t) G Baya(s,t), Yo(s, t) = Biy\(s, t), 

y3(s, t) = Bay4(s, t), ya(s, t) = B3y3(s, t) + Biyi(s, t), (2 1) 
yi(1, t) = —aıy2(1,t), y2(0,t) =0, l 
y3(1, t) = —agya(1,t), ya(0, t) = 0, 


in which 6; = K, 62 = p~', 63 = EI and 64 = I>". Here and Hereafter, the dot and the prime 
represent derivatives with respect to time and spacial variables, respectively. Set Y(s,t) = (yi(s,t), 
yo(s,t), y3(s,t),y4(s,t))' and the first four equations of (2.1) can be written as port-Hamiltonian 
system 

Y(s,t) = P\[HY(s, t)]' + PoHY(s, t), (2.2) 


in which H = diag((1, Bo, p3, Ba) and 


000 -1 010 0 
000 0 100 0 
=o 9 0 of 2 ~|0 0 04 
100 0 001 0 


The state space of (2.1) or (2.2) is X := L?((0,1];C*) with inner product 
1 
(Y, Z}x = f Y(s)*HZ(s)ds, V Y,Z €X. (2.3) 
0 


The rest of this subsection is employed from [3]. To formulate the boundary conditions of (2.1) or 
(2.2) in a better manner, we introduce the boundary effort and boundary flow, which are defined as 


1 
eð = — 
eye) 


respectively. Hence, we consider the operator 


H[Y (1) + ¥(0)] and fo = BPHYA) ~¥(0)], (2.4) 


AY (s) = P,(HY(s))’ + Po(HY(s)), Y Y(s) € D(A) (2.5) 


on the state space X and the domain 


D(A) = [Y EX: HY € M (0,1; We £] =o}, (2.6) 


in which H+([0, 1]; C4) is Sobolev space of one order and Wg € C4*® equals 


a, BB} 0 0 BB a 0 0 

ij- 0 0 0 0 

v2 | 0 0 az bab3’ 0 0 Bibs! az 
0 0 —1 0 0 0 0 1 


Wp = 


Thus, the abstract differential equation 
Y(s,t) = AY(s,t), Y(s,0)€ X (2.7) 


is equivalent to (2.1). To obtain the main result of this subsection, we introduce two lemmas. 
0 Ll 


Lemma 2.1 Let J, be the identity operator on C4 and © = i 0 
4 


| then the matrix Wg has 


full rank and satisfies 


WpmuW3 > 0. (2.8) 
Proof. Let Wg = [W1, W2] and 
a, B2B} 0 0 B267! ay 0 0 
1 —1 0 0 0 1 0 1 0 0 
W, = — = W: = —= — , 
1 AZ| O0 0 aœ Bp! °? A| 0 0 BB œ 
0 0 —1 0 0 (0) 0 1 
then we have 
WB EWS = WiWs + WW 
2018871 —BoBy" 0 0 2018287! 267" 0 0 
_ 1] B67" 0 0 0 ii — b267 ' 0 0 0 
2 0 0 2028483" -bab3 | ° 2 0 0 2028483' bab3" 
0 0 BaPs * 0 0 0 — BBs" 0 


= diag(2a1 b287", 0, 2028483", 0) = 0, 


which implies (2.8) holds. 


The following lemma comes from [6, Lemma 7.2.1]. 
Lemma 2.2 For any Y € D(A), the following results hold 


Re (AY, Yx = —a bı b2ly2(1, t)? = a23384\ya(1, t)’. (2.9) 


Now we can show that (2.7), and hence (2.1), is exponential stable. More precisely, 


Theorem 2.1 The operator A generates a Co-semigroup of contractions T(t) and T(t) is expo- 
nential stable. 


Proof. Theorem 7.2.4 of [6] and Lemma 2.1 imply that the first statement is right and Theorem 9.1.3 
of [6] and Lemma 2.2 mean that the second statement is right. 


2.2. Uniform exponential stability of (2.1) 


In this section we discuss the uniform exponential stability of Timoshenko beam (1.1)-(1.3). For 
this purpose, we perform the spacial semi-discretization for (2.1). Let N € N be an integer and 
h =1/(N +1) be step size, insert N + 2 points and N + 1 points, which are denoted by s; = ih(i = 
0,1,---,N+1) and v; = (j +1/2) h(j = 0,1,--- , N), respectively, in the domain 0, 1]. f; only be 
the value of a function f(s) at node s; = ih(i = 0,1,---,N+ 1). This means that y, ;(t) denotes 
yi(si,t) for l = 1,2,3,4. The notations ôs fj}; = (fj41 — fj)/h and fj,1 = (fj41 + fj)/2 denote the 
central divided difference operator of first-order derivative at f,(v;) and the average operator at f(v,), 
respectively. 

Now we are in a position to give the spacial semi-discretizaiton scheme of (2.1). The independent 
variable s in equations (2.1) is replace by vj, i.e., 


oe yo(v;,t) = biyi (vz, t), 
ys (vj, t) = Baya (vj, t), Ya(vj,t) = Bsy3(vj,t) + Bryr(v;,t). 


Approximating the spacial derivative by difference operator and applying the average operator to the 
time derivative respectively, we get 


vı +h 1 (t) = BodsYo jth a{t)— Baya j+ t), 
Yoj+h i(t) = b18 sYı, j+s i(t), J =0,1,-: N, 
Y3,j+3 (t) = BadsYa +h 1(t), (2.10) 
Vaji (t) = P3bs¥3,j42(t) + Pita j+ (Ë), 
yı, n+ (t) = onya Waal) y2,0(t) = 0, 
y3,n+i(t) = —a2ya,n+i(t), y4,o(t) = 0, 
by considering the boundary conditions in (2.1). 
In the state space C4ON+1) define 
Z 41,1 Zk,0 
Za, 21,2 Žk,1 Žlj»Žkj E€ C, 
Xr = 4 Zr = z e CAN) : Zin = ; , Zehr = f ; 
oo k =1,3,l = 2,4, 
4h Z1 N+1 Žk,N 
with the inner product 
= i N 7 Z 
(Zh, Zh)h = hiai Žj=0 (Bizitu) ,VZh, Zh € X}, (2.11) 
22,0 = 24,0 = 0, Z1, N+1 = —Q122,N41, 23,N+1 = —Q224,N4+1- 


The state variable of (2.10) is 


Yn (t) = (Yin(t), Yon (t), Yan(t), Yan (©) 


with 
Yin(t) = (y1, (t) yin (t)), Yaon(t) = (Y2,1(4),-++ y2,n4i(t)), 
Y3n(t) = (y3,0 (t) Y3, N (t)), Yan(t) = (ya (t) Ya, N41 (t) 


then the energy En (t) is 


1 » oN 
En(t) = z nC), Yn ))x, = 5 228 ACh ged i (2.12) 
1=1 j=0 


We have a discrete counterpart of Lemma 2.2 


Lemma 2.3 The solution Y, (t) to the system (2.10) satisfies the balance equations of discrete 
version: 


En(t) = -a1 b1 boly, N41 (H|? — a2b3balya, n1 (t). (2.13) 


Proof. Multiplying the first four differential equations of (2.10) by hSiy;;41(t) (l = 1,2,3,4) and 
adding them together, we obtain 


N 
E, (Y (Y; t) =hb1b2 X uijt )OsYo,j42 1(t) + 


=0 


+h (t)ðsY1 j1 (t) 


2 


rs 


N 
+hp3p4 Yo ujal )osYaj4a(t) + 
j=0 j=0 


Me IM 
S S 
rs a 


4 (t)ôsY3 j+ (t) 


However, a simple calculation shows that 


i So nseat )5syo,544 ( D+ Dorsal Jõsy1 j4 (t) | = y1, N+ (t)y2, N41 (t) — y1,0(t)y2,0(t), 


N 


h 5 Y3,j44 (t)ds¥a,j44 (t) + 2 Ya j+ (t)ôsYs jti (E) | = Ys, w+ (t)ya, n+ (t) — y3,0(t)ya,o(t). 
j=0 


The above two identities and boundary conditions imply that the equality (2.13) holds. 


Moreover, we additionally need the following lemma (see [20-22]) to obtain main result of this 
subsection. 


Lemma 2.4 Let {uj} 4", {u; 44) and {w:}}" be sequences consisting of real numbers, then 
we have 
N h2 
h J CEN + Uj 410s; 41 Wip + Uj Upg is Wip + 7 Oslin £O8Vi 41 95Wi4 2 
i=0 


=UN+1UN+1WN+41 — UNUVOWO- 


To represent (2.10) into vectorial form, we introduce some matrices. Let An, Bh, Ch, and Dp be 
square matrices of order N + 1: 


'202301.00192v1 


chinaXiv 


Ch = diag(0,--- ,0,-a1), Dn = diag(0,--- ,0,—az2), respectively. Set 


Hn = diag(61Iw41, b2IN+1;, B3In41, Bain +1), 
Ap, 2-'C), 0 0 


Pa = 0 0 An 271 Dy i 
0 0 0 A} 
0 b2B, 0 —ß4 AF 
y, = —2, Bn 0 0 0 
E 0 0 0 BaB J’ 
Bi An 0 —B3 Bp 0 
0 0 0 0 
Q = 0 Bh-1C, 0 0 
h= 0 0 0 0 
0 2-18,Ch 0 B3h-1Dp 
After these preparations, we can rewrite (2.10) as 
Yn (t) = AnYn(t), Yn (0) € Xn, (2.14) 


where An = ©, (Yn + Qh) because ©), is evidently invertible. 

It is easy to see that ®, is corresponding to the average operator of time derivative of (2.10). If 
one replace ® by identity operator, then the classical finite difference scheme of (2.1) is easily restored 
from (2.14), i.e., 

Ya(t) = (Wn + On)¥n(t), Ya(0) € Xn. (2.15) 


Here we explain the significance of the discrete scheme (2.14). We plot two figures in Figures 1 and 
2, respectively. Figure 1 depicts the maximal real parts of the eigenvalues of our discrete scheme (2.14) 
and the classical semi-discrete scheme (2.15) for N = 10:5: 400. Figure 2 depicts the distributions of 
the eigenvalues of (2.14) and (2.15) in which N = 100. We see that the real parts of the eigenvalues 


T 
x without the average operator x without the average operator 
gar x with the average operator | 400 F x with the average operator 


-100 F J 
-200 + 


-300 F 


Maximal real parts of eigenvalues 
i : pi i : 
a 
Imaginary part 


-400 F 


-0.8 -500 


0 50 100 150 200 250 300 350 400 3 25 2 15 ali -0.5 0 


N Real part 
(a) Fig.1. Maximal real parts of eigenvalues (b) Fig.2. Eigenvalue distributions 


of (2.15) approach to zero and those of (2.14) approach to a negative number from two figures. In 
both figures, we take 8; = b2 = 63 = b4 = 1, ay = ag = 1. Numerical simulation results show that 
(2.15) is not uniformly exponentially stable. This conclusion is consistent with earlier researches of 
[16]. However, Figure 1 manifests (2.14) is perhaps uniformly exponentially stable. In the remaining 
part of this subcection, we give strict proof about uniform exponential stability of (2.14). 


As announced in the Introduction, we will apply the perturbation theory of semigroup to discuss 
the uniform exponential stability of (2.10) or (2.14). So we first study the uniform exponential stability 
of the system 


Ui j+} (t) = P2ðsY2 j+4 (t), U2 j+} (t) = Biðsyı jt i(t), 
U3 j+} (t) = BadsYa,j44 (t), Va j+ (t) = B30sY3, „j+ a(t J (2.16) 
yi Nil) = —arye,nei(t), ys,N+1(t) = —azya N+1(¢), 
y2,o(t) = y4,0(t) = 0, j=0,1,:-- N. 
we also represent (2.16) into vectorial form, let 
0 BoB} 0 0 
ï, = — bı Bp, 0 0 0 
iT 0 0 0 BaB} 
0 0 —b3Br 0 
Then we can write (2.16) as 
Yn (t) = AnYn(t), Y;,(0) E Xa}, (2.17) 


where An = 3 (Yr +9). Set Mı = max{1/V 6162, 1/V B361} and the uniform exponential stability 
of (2.16) or (2.17) is presented as follows. 

Theorem 2.2 For any initial value Y, (0) € X}, the energy Ep (t) of the solution to (2.16) decays 
uniformly exponentially with respect to the step size, i.e., there exist two positive constants M and 7 
independent of h such that 


l+eM, E 


En(t < Me ™*E ith M = ? = 
a ) _ E n (0), wi 1— eM “ l+eM,’ 


(2.18) 


for0<¢e<1/Mj. 


Proof. Firstly, we introduce auxiliary functional 


N 


galt) =AY eritir eases © + 54200543 ana À 
j=0 


and Lyapunov functional G;,(t) = En(t) + evn(t). It is easy to see that |y;,(t)| < M1 En (t) and hence 
one has 
(1 — Mi) Ep (t) < Galt) < (1 + €M1)En(t), (2.19) 


in which the condition 0 < € < 1/M; ensures that the Lyapunov functional G, (t) is positive definite. 
Secondly, differentiating y,(t) with respect to t and applying (2.16), we get 


N 


pnt) =h 7 [Prsni sY j+ (t)Y, j3 Ct) + B28j445s¥29+3 Ewa (8) 
j= 


N 
hX [8355.4 105Y3.54 1(t)Ys. 54 1(t) oF ETE EOE] f (2.20) 
j=0 


Applying Lemma 2.4 to each term of the right-hand side of (2.20), we have 


4 


ok 
Gn(t) = —En(t SDT j O)? 
l=1 j=0 
= (b2 F 6,04) lyo,n4i(t)|? + = 5 (pa F 305) lya, n4ilt)|?. (2.21) 


At last, differentiating the Lyapunov functional G, (t) and applying Lemma 2.3, we have 
Gn (t) = Én (t) + ep, (t) < —nG; (t), 


since € < 2a1/162/(B1a7 + b2) < VBiB2 and e < 2@2ß364/(Ba + a383) < V Bsba. Thus, applying the 
comparison principle (see [22] or [31, Section 2.3]) and (2.19), we obtain (2.18). 


Furthermore, we define the disturbance operators Dp as 
Dn Yn (t) = (—B4An Yan (t), 0, 0, 81 AnYin(t))" 


Then we have the following theorem. 
Theorem 2.3 Ife is the same as in the last theorem and 2V M 184 < 7 holds, then the system 
(2.10) or (2.14) is uniformly exponentially stable in the sense of Theorem 2.2. 


Proof. Denote the Co-semigroups generated by the operators An and Ap as Th (t) and Ton (t), respec- 
tively. We should know the norms of the operators D, to analyze the exponential stability of the 
semigroups T},(t) by the theory of bounded perturbations of Co-semigroups [6, Theorem 10.3.1]. In 
fact, by the definitions of the operator Dp , it is easy to see that for any Zp € Xp, 


N 
Pn Zalk, =A [Balaza 543 (OI + BalBr2rg4 3)? 
j=0 


It is not difficult to verify that ||D),||%, = 8184 holds. From Theorem 2.2, we know Ton (t) are uniformly 
exponentially stable since 2E;,(Y;t) = ||Tn(t)Y;,(0)||?, for any t > 0. Whereas (2.18) implies that 


I|Ton (t) Yn (0)| 


xX, = 2En(t) < 2Me~™ Ep,(0) = Me~™||¥,(0) ||, Y Yn(0) € Xn, 
which is equivalent to ||Tpn(t)||x, < VMe™™/2. Using [6, Theorem 10.3.1], we obtain 
ITa Oln, < VM VMP), 


The condition 2y M8164 < n ensures —n/2+ V M||Dh]|x, < 0 which implies that the semigroups Ty (t) 
are uniformly exponentially stable. 


Remark 2.1 One may wonder whether or not there is € such that all the conditions given in 
Theorem 2.2 and Theorem 2.3 are satisfied. In fact, let 6182 = 8384 = 1 hold, then one should choose 
6, and £4 to satisfy 

é"(1—e) 


A(1 +e)? 


which ensure that 0 < € < 1/Mı and 2V M8164 < 7 hold. This manifests that 6;(1 = 1,2,3,4) is not 
always suitable for the conditions given in Theorem 2.2 and Theorem 2.3, this is a restriction of this 
method. 

Remark 2.2 Theorem 9.1.3 of [6], which is used to prove Theorem 2.1, is very effective to study 
the exponential stability of port-Hamiltonian system on infinite-dimensional space. To the best of our 
knowledge, there is no counterpart for uniform exponential stability of a family of port-Hamiltonian 
systems on finite-dimensional spaces. From the proof of Theorem 9.1.3 of [6](see also [3]), we know that 
the exponential stability highly relies on a kind of final state observability of a dissipative boundary 
control system. Thus, how to generalize this result to testify uniform final state observability of a 
family of dissipative boundary control systems is an open problem. 


BiBa < with O<e<1, 


3. Related results of the beam (1.4) 


In this section, we study the uniform exponential stability of Timoshenko beam (1.4) with interior 
damping. Two kinds of methods applied in this section are completely different in contrast with the 
methods of last section. 


3.1. Exponential stability of continuous beam (1.4) 


Using the auxiliary functions y(s,t)(1 = 1,2,3,4) of subsection 2.1, we can transform (1.4) into 
the following system 


yi(s,t) = BoyS(s,t) — Baya(s, t), 
Yyo(s,t) = Pry (s, t) — B2y2(8,t), 
¥3(s,t) = Baya (s, t), (3.1) 
jals, t) = Bsy3(s,t) + Piyi(s,t) — Baya(s, t), 
1(1, t) = yo(0, t) = ys(1, t) = y4(0,t) = 0 
The first four equations of (3.1) can be rewritten as the form 
Y (s,t) = Pi[HY(s,t)]’ + [Po + PIHY(s, t), (3.2) 
in which H, Po, and P, are defined in subsection 2.1 and P = diag(0, —1,0,—1). 
Now, we consider the operator 
AoY (s) = Pi(HY(s))’ + PHY (s), V Y (s) € D(Ao) (3.3) 
on the state space X and the domain 
D(Ao) = fy €X: HY € H ([0, 1]; C$), W 3 || = o} , (3.4) 


in which W8 € C4*8 is derived by setting a; and ag to be zero in Wp. 
Let the operator B be defined by BY (s) = PHY (s), for all Y (s) € X, then the abstract differential 
equation 
Y(s,t) = (Ao + B)Y(s,t), Y(s,0) €X (3.5) 


is equivalent to (3.1). The operators Ao, B, and Ao + B satisfy the following properties. 
Property 3.1 For any Y € D(Ao), we have 


Re( oY, Y) = -f B3lyo(s,t)|? + 22lya (s, t)|?ds. (3.6) 


Proof. In fact, for Y € D(Ap), it follows from the definition of Ag and the inner product (-,-)y given 
in last subsection that 


1 
(AoY, Y)x =f Bıb2y5(8, t91 (8, t) + b1b291 (8, t928, t) — B3|y2(s, t)|?ds 
1 
+ f B3 Bay4(s, t)¥3(s,t) + B3 bays (s, t)y4(s, t) — Bi lyals, t)|?ds 
and 
1 pen —* 
(Y, AoY)x =f B1b24'2(5, t)y1 (s, t) + Bi B2y'1(s, t)y2ls, t) — B3|yo(s, t)|?ds 


1 
+f B3Bay' 4(s, t)ys(s, t) + B384y's(s, t)ya(s, t) — B3|yo(s, t)|?ds. 
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Adding two identities above together, we derive 
Re(AoY, Y)x Ti Baboly2(s, t)Y1 (s, t) + y1 (s, t)¥o(s, t) tas- f b3ly2(s, t)? + Balya(s, t)|?ds 


+5 [ Padalvals.105(6st) + vols, Dals, ds, 


which yields (3.6) by W$ [f3, e5] =0 & y1(1, t) = y2(0, t) = ys(1, t) = ya(0, t) = 0. 


Property 3.2 The operator B is skew-adjoint and hence Re(AoY,Y)x = Re((Ao + B)Y,Y)x < 
0. 


Proof. We only prove the first statement since Re(BY, Y)x = 0 if operator B is skew-adjoint. Whereas, 
it is easy to see that 


(BY, Z)x + (Y, BZ)x 
1 


1 
=P p4 f —y4(s, ty, (s, t) + YL (s, t)y4(s, t)ds F Biba f —y1(s, t)y4(s, t) + ya(s, t)yı (s, t)ds 
=0, 


which shows that B* = — B. 


Properties 3.1 and 3.2 imply that the operators Ao and Ag+ B are both dissipative and they further 
imply that the point spectral sets 0,(Ao) and op(Ao + B) of Ag and Ao + B are contained in the closed 
left half-plane of C. Moreover, we have more stronger results. 

Property 3.3 The point spectral sets o,(Ao) and o,(Ap + B) of Ao and Ap + B are contained 
in the open left half-plane of C. 


Proof. We only prove o,(Ao) is contained in the open left half-plane of C since another proof is similar. 
To this end, we show now by a contradiction argument. If there exist 6 € R and nonzero Y € X such 
that AoY = ißY , then it follows from (3.3) and (3.4) that 


iByi(s) = Boyo(s), ißy2(s) = Biyi (s) — Baye(s), 
iBy3(s) = Bay4(s), *Bya(s) = Bsy3(s) — Baya(s), (3.7) 
yi(1) = y2(0) =0, ys(1) = ya(0) = 0. 

It follows from Re(AgY, Y)x = Re(i8Y, Y)x = 0 and Property 3.1 that yo(s) = y4(s) = 0. Thus, (3.7) 


is reduced to 
{ iBy1(s) = y1 (s) = tBys(s) = y3(s) = 0, 
yi(1) = yo(0) = ys(1) = ya(0) = 0. 
Poincaré inequality, y,(1) = y3(1) = 0 and yj (s) = y3(s) = 0 imply that yı(s) = y3(s) = 0. This leads 
to a contradiction and thus o,(Ao) N iR = Ø holds. 


Property 3.4 Let f(s) = 6;'5'(s? + b28) and g(s) = 6;'8;,'(s? + 64s) be polynomials and 
Un = 7/2+ nn (n € Z), then the point spectral set op(Ao) is given by 


op(Ao) = {A: FA) = =u n EZ} U {A : gA) = =u n EZ}, (3.8) 


and the corresponding eigenvectors form an orthonormal basis. 
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Proof. Assume Y € D(Ap) such that AoY = AY, then it follows from (3.3) and (3.4) that 


Ayi(s) = Boyo(s), Aya(s) = Fiyi(s) — Baye(s), 
Ays(s) = Baya(s), Aya(s) = Bsy3(s) — Baya(s), 
yi(1) = y2(0) = 0, y3(1) = y4(0) = 0. 


The above equations are then equivalent to 


yo (s) = fA)y2(s), ya (s) = g(A)ya(s), 
Boy(s)=Ayi(s),  Bay4(s) = Ays(s), 
yi(1) = y2(0)=0,  y3(1) = ya(0) = 0. 


yo(s) and y4(s) can be easily solved by y2(0) = 0 and y4(0) = 0: 


yo(s) = Cı sinh(y f(A)s), ya(s) = C2 sinh(vg()s), 


where Cı and C2 are a sae determined later. Moreover, yı(1) = 0 and y3(1) = 0 implies 
fA) = iun and \/g(A) = iun (n € Z). Thus, this means that (3.8) holds. More precisely, we obtain 
the point spectral set a consisting of two branches o (Ao) and a? (Ao), which are respectively 


given by 
o}(Ap) = fa: _ z+ la 481 Bop? | , 02(Ap) = fa _ zb vii- 483 ban, | ned. 


The eigenvectors Y; and Z> corresponding to eigenvalues A% and y% are 


n 


4 1 1 
Y7 = Cin - c0s(tns) sins}, 0,0 and Z> = Con, fo. 0,-—= cos(Hns) sin) : 
n Yn 


n 


respectively. The parameters Cin and Con are chosen such that ||Y,*||x = 1 and ||Z*||x = 1. It follows 
from basic theory of function that {Y= : n € Z} U {ZF : n € Z} forms orthonormal basis of X. 


Under these preparations, we can prove an exponential stability result on the operator Ag. 
Theorem 3.1 The operator Ap generates a contractive semigroup which is exponentially stable. 


Proof. We first show that Ag generates a contractive semigroup. By Property 3.1, we know that Ag is 
dissipative. Thus, it is sufficient to to show that Ran(J — Ao) = X by Lumer-Phillips Theorem. That 
is to say, for each Z € X we have to find Y € D(Ao) such that (I — Ao) Y = Z, which is equivalent to 


yi(s) — Boya(s) = 21(s), (1+ B2)yo(s) — Bry (s) = 22(s), 
ys(s) — Baya(s) = z3(s), (1+ Ba)ya(s) — Bays(s) = 24(s), 


By the same method of [6, Example 6.1.10], we can obtain unique Y given by 


( As ) 7 ( ba ees ) + f B(s- r)(a(s = 7), 22(s- 7))"dr, 
( ) ( (ofan iia ahe) ) 7 [ $2(s—r)(2a(s— 7), za(s = 7))' dr, 


in which y1 = 4/87 83 (1 + b2), Y2 = 1/63 BI (1 + Ba); 
1 1 
--a | sinh(y1(s — T))z1(s — T) + Boy sinh(y1(s — T))z2(s — 7)dr, 
0 
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69 = —— : T f ion cosh(72(s — T))z3(s — T) + cosh(y2(s — T))z4(s — T)dr, 


sinh( 
(s) = sinh(75) 6271 sinh(718) 
g (28o71)~1 cosh(718) cosh(y1s) J’ 


= sinh(72s) Baya sinh(728) 
=) — Cre cosh(728) er ) 


On the one hand, this means that (I — Ao) is bijective and hence Ao generates a contractive semigroup. 
On the other hand, it follows from [32, Section 22.2] that (I — Ao)~! is a compact operator. Therefore, 
the spectral set o(Ao) and the point spectral set o,(Ao) coincide and Property 3.4 further implies 
spectrum determined growth assumption holds. The exponential stability of the semi-group generated 
by Ap is finally derived since Property 3.4 manifests that the spectral bound of Ag is negative. 


Applying this theorem and [13, Theorem 5.1], we can prove the exponential stability of the semi- 
group generated by the operator Ag + B. To this end, we introduce the definition of B is A-compact 
since it is involved in [13, Theorem 5.1]. 

Definition 3.1 Let A and B defined on Banach space X be two linear operators such that 
D(A) c D(B), we shall say that B is A-compact, if for any sequence (£n)n C D(A) with both (£n)n 
and (Az,,), bounded, (Bx,), contains a convergent subsequence. 

Theorem 3.2 The operator Ag + B generates an exponentially stable semigroup. 


Proof. Theorem 5.1 of [13] asserts that the semigroup generated by Ag + B is exponential stable if and 
only if op(Ao) lies in the open left-half plane if B is Ag-compact. But, we have proven that o,(Ao) is 
contained in the open left-half plane in Property 3.3. It is sufficient to verify that B is Ag-compact. 
Thus, let (Yn (s, t))n C D(A) with both (Yn(s, t))n and (ApYn(s,t))n bounded, we have (yin(s,t))n and 
(Yin(S,t))n are uniformly bounded for l = 1,2,3,4 by the definition of the operator Ag. Combining 
Theorem 2 of [32, Section 22.1], which is Rellich’s standard of precompact set in L?(Q) with Q being 
the open set of R”, and BY = (—B4y4(s, t), 0,0, B1yi(s,t)), we know that B is Ag-compact. 


3.2. Uniform exponential stability of (3.1) 


To study the uniform exponential stability of (3.1), we discretize it by the average central-difference 
method of section 2.2. It is easy to know that the semi-discretized system of (3.1) is 


Ùi j4 i(t)= B2ðz W jth i(t) — Bawa j+ (t), t>0, 

Wo j4 i(t)= Bi dx Wy ge (t) — Bowe 542 (t), j =0,1,--+ ,N, 

3 541 (t) = Badew4 542 (t), (3.9) 
Wa j+ (t) = B3d2W3,542(t) + biw 541 (t) — Bawa ja (t), 

wi N41(t) = wa,o(t) = ws,n41(t) = wao(t) = 0. 


Let Wn (t) = (Win, Wan, W3n, Wan)! be unknown variables consisting of 
Win(t) = (wio(t), +, wiw(t)), Wan(t) = (war(t), =, w2,N41(t)), 


Wsn(t) = (wao(t), =+, ws,w(t)), Wan(t) = (war(t), =, waw+i(t)). 


An, Bn, and Hn are given in the last section, we therefore define the matrix 


0 B} 0 —A} 

o| -B, -A4 0 0 
Gr=| o 0 0 B’ 

A, 0 -Br —Aj 
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Using these matrices, we rewrite (3.9) into a vector form 
®,W,,(t) = GrH Wn (t). (3.10) 


The state spaces of (3.9) is still denoted by X, which is the linear space C4(N+) equipped with the 
inner product (Wh, Va), = h(®nWnh, Hn ®nVn) for Wi, Vn € Xn, in which (-,-) is natural inner product 
of C4N+) or CNHI, The classical finite difference scheme of (3.1) is easily obtained from (3.10), 


Wn (t) = GhHnrWa (t). (3.11) 


Similarly as in subsection 2.2, we also plot two figures to explain the effective of our numerical 
approximating scheme (3.10). Figure 3 depicts the maximal real parts of the eigenvalues of our discrete 
scheme (3.10) and the classical semi-discrete scheme (3.11) for N = 10 : 5 : 400. Figure 4 depicts the 
distributions of the eigenvalues of (3.10) and (3.11) in which N = 100. In both figures, we take 
6, = 10, b2 = 20, B3 = 30, 64 = 40. We can draw the same conclusion with those of subsection 2.2. 


That is to say (3.10) is uniformly exponentially stable and (3.11) does not possess this property. 
q x10 
x without the average operator x without the average operator 
x with the average operator 0.8 F x with the average operator 
g Pf 0.6 
2 
oO 
2 al 0.4 H % 
2 e 
oO I 4 
z & 0.2 7 
£ -6 z | a eee 
2 > x 
è E -0.2% | 
g S | 0.4} * 
E] 
= 0.6 F 
-10F 
0.8 F 
-12 1 f i 1 1 1 f 4 1 1 f 
0 50 100 150 200 250 300 350 400 -20 -15 -10 5 0 
N Real part 
(c) Fig.3. Maximal real parts of eigenvalues (d) Fig.4. Eigenvalue distributions 


To prove (3.10) is uniformly exponentially stable, we firstly show that Ay = Dr GhHn is dissipative 
for every step size h. 

Lemma 3.1 The operator Ap := ®,'GrHn is dissipative on the space Xj, i.e., for h € (0,1) and 
Yn E€ Xp, 


Re(AnYn, Yn)x, = —O2h(A] Yon, An Yan) — BZh(Ap Yan, Ap Yan) < 0. (3.12) 
Proof. Here and hereafter, we assume Yp = (Yin, Yon, ¥3n, Yan)! € Xp with 


Yin(t) = (y1,0 (t) , yin (t)), Yont) = (y2,1 (t), , Y2,N-Ṣ 
Y3n(t) = (y3,0 (t) Y3, N (t)), Yan (t) = (ya, (t), Y4, NṢ 


T T 

=. e 
—_~— 
bho 
Wn 


but we need to assign some imaginary components such as yi, +1(t) = y2,0(t) = y3,w41(t) = y4, (t) = 
0. By the definition of the operator A, and the definition of the inner product (-,-)x,, we have 
(An¥n, Ya)x, = h(GnHiYn, Hnn Yn) 
=61B82h(By Yon, AnYin) — 6182h(BuYin, An Yan) — B2h(Ap Yon, An Yon) 
+ B3Bah{ By Yan, AnY3n) — 838ah(BnYan, An Yan) — Bth(Ap Yan, An Yan)- (3.13) 
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and 


(Yh, AnYn)x, = h(GnYn, HhGh Hn Yn) 
=B1Boh(AnYin, B} Yon) — b1b2hl A} Yon, BhYin) — BZh( Aj; Yon, Ap Yon) 
+ B3Bahl AhYsn, B} Yan) — B3Bah( AF Yan, BhYsn) — BARLA Yan, Ap Yan). (3.14) 


By some simple calculations, we know that 


h( By Yon, AnYin) — h( Aj, Yon, BhYin) 


N = _ N = _ 
=> (yo,j+1 — Y2,5) (Vi joa + 91,5) = (yi j+1 — Yi) (Goji + Y2) 

; 2 4 2 

j=0 j=0 
=Y1,N+192,N41 — Y1,092,0 = 0. (3.15) 


Similarly, we also have 


—h(BnYin, Ap, Yon) + h(AnYin, Bp Yon) = 0 
h(BY Yan, AnYan) — hA; Yan, BaYan) = 0. (3.16) 


—h(BnY3n, Ay, Yan) + h(AnY3n, By, Yan) = 0 


Thus, (3.13)-(3.16) implies that (3.12) is right since (AnYp,Yn)x, + (Yn, AnYn)x, = 2Re(AnYn, Yn)x,- 


From last lemma, we know that An generates a family of semigroups of contraction Ta (t). In 
the sequel, we apply frequency standard of uniform exponential stability for a family of contractive 
semigroups, which is given in [15], [29] or [30], to study the uniform exponential stability of Tp (t) for 
all h € (0,1). 

Theorem 3.3 Let h* > 0 and (.S;,(t)) be a family of semigroups of contraction on the Hilbert 
space (Xn), and let (An) be the corresponding infinitesimal generators. The family (.S;,(¢)) is uniformly 
exponentially stable if and only if the two following conditions are satisfied: 

(i) For all h € (0,h*), iR is contained in the resolvent set p(Ap) of (Ap). 

(ii) suPpe(o,nx),ger IGBI — Anlaa < ©. 

We first show that (i) holds for the operator Ap. 

Lemma 3.2 For all h € (0,1), {48, B € R} C p(An). 


Proof. Suppose this conclusion is false, i.e., iĝ € o(An), then there exists 0 Æ Yp € Xn and 8 € R such 
that AnYn = iBY),. On the one hand, it is easy to see that Re(AnYn, Yn)x, = Re(iBYn, Yn)x, = 0. By 
(3.12), we know that (Aj Yan, A}, Yon) = 0 and (Aj Yan, A} Yan) = 0, which implies that Yo, = Yan = 0, 
since Bp is invertible. On the other hand, ALY; = ißYn is reduced to 


tBAnYin =0, — b1BhYin = 0, Bi: AnYin — P3BnY3n = 0, 


which yield that Yi, = Ysa = 0 since By, are invertible. Thus Y, = 0 is derived and this contradicts 
Yn Æ 0. 


Now we can present the main result of this subsection. 
Theorem 3.4 The family (T, (t)) is uniformly exponentially stable, that is to say there exist two 
constants M > 0 and w > 0 (independent of h € (0,1)) such that 


ITO) < Me™, Vt > 0. (3.17) 
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Proof The proof is based on Theorem 3.3. Notice first that, for all h € (0,1), the family (T),(¢)) 
form a family of contraction semigroups (see Lemma 3.1). The fact that the family (Aj) satisfies 
condition (i) follows from Lemma 3.2. In order to show that the family (Ap) satisfies condition (ii) we 
use a contradiction argument. If the condition (ii) is false, then there exist b, € R, hn € (0,1), and 
Zh, © Xh, Such that 


I|(@bnTh, — Ån, ZR, IIx, > runba NZR, Il, 


in which vp, = |ibn +max{ b2, S4}|?, n = 1,2,---, and the norm ||-||x, is induced by the inner product 
of Xp,,. Let (tbp Jn, — An.) ZR, = Fy’ and we have 


IER xn, Z NVnhp || libnInn, — Ang) ER, [kin 
Setting YP = BP [ER Xn we obtain Ye xa, = ] and 
[UR IIx,,, <n tvs hn, with UR, = (ibnIn, — Ån, Y; . (3.18) 
More precisely, for Ug, = (Uf), , Udh,, o U3h, Uth, ) we have 


OR lkn, SAn (61) Arn Oth, I? + Bell An,, Ozh, I? + B3llAn, Unn I? + BallAn,, Oth, l] 
<u n he. (3.19) 


It follows from Cauchy-Schwartz inequality, (3.12), and (3.18) that 


Re (UR, Yi Jx, = -Re (An Yi YR Y 


hn 
=hn [83| An, Yohn N? + Ball An, Yangll?] < vn nhn. (3.20) 
This means that 
hnll An, Yan (I? = OW, hin); hnll A}, Yan, lI? = O(v nhn), (3.21) 


in which the notation fa = O(n~+) denotes that there is a positive constant C (independent of n) such 
that ||fn]| < Cn7* for n € Z+. Now we consider the identity Ug, = (ibnIn,, — An, ) Y}, which can be 
expanded as 


ibn AnY ih, — G2By Yah, + B4An Yah, =AhUih, 
Pi BrY th, + (ibn + B2) Ap Yoh, =An Uh, 
ibn AnY3h,, — PaBn Yah, =AnU3h,,, 


— Bi AnYin, — B3 BhY3h, + (ibn + Ba)An Yah, =An Uih, 


n 


More precisely, we have 


ibay jga — B2dn¥> 543 + Badr yy ; 1 =u]; 1; (3.22) 

(ibn + B2)Y3 j4} — Prde¥y jti SUZ 5425 (3.23) 

ibnY3 543 — Bade Vy j I =u; j 1 (3.24) 

(ibn + Bayt 544 = B3ba¥3 542 — Biyy jpa H =01 Ug 544, (3.25) 


for j =0,1,--- , Nn and yf y 41 = Yeo = YEN, 41 = Yio = 9. It follows from (3.23) that 


AYN, = hnus N 43 — hn (ibn + b2)Y3 Nath? 
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holds by setting j = Nn. (3.19) and (3.21) further imply that 


Win, < lt + Pahl wy, gal? + ee 2 < Clidy + max{ 32, ba}l?vz n thn 
which can be expressed as 
lvi Na l” = O(n Rn), 
since Vp, = libn + max{ 62, B4}|?. Similarly, by setting j = Nn — 1 in (3.23), we get 
Biytn,-1 = Pi¥in, — Rn(ibn + B2)Y3 N, -3 + And 2U N,- 
and 
lvi Nail = O(n *hn). 
We can repeat this process to obtain 
luis = O(n™'hn), for all j = 0,1,- , Nn, 


which implies that 
Nn 
hall AnYih, I? = hn Do Ivy pal? = O(n), (3.26) 
j=0 


since hn(Nn +1) = 1. Using (3.25) and the same method as above, we know that Y3}, also satisfies 
lyzl? = O(n ha), for all j =0,1,--- , Nn, 


and 


Nn 
hnll An Yh, l? = hn X lv jag)? = O(n"). (3.27) 
j=0 


Finally, it follows from that (3.21), (3.26), and (3.27) that 


Yi Ik, =hn [81ll An, Ying ll? + Ball Ah, Yon, II? + Bll Ann Ysha I? + Ball An, Yah, 


*] =O), 


which leads to a contradiction since ||Y;? ||x,, = 1 hold for all positive integers n. The proof of the 
Theorem 3.4 is complete. 


4. Concluding remarks 


Timoshenko beam is a basic vibration model and plays important role in engineering. This paper is 
devoted to uniformly exponentially stable approximations for the Timoshenko beam with both interior 
damping and boundary damping. This means that we study it from the viewpoints of the numerical 
approximating and control theory. It is well-known that there are many discretization methods to 
discretize the spatial variables. It is not an easy job to pick one which preserves exponential stability 
among so many semi-discretization methods. Moreover, Timoshenko beam is described by two coupled 
wave equation and the crossed terms bring troubles in verifying the uniform exponential stability. To 
bypass these challenges, perturbation of exponentially stable Co-semigroup is introduced and uniform 
exponential stability of both continuous and semi-discretization systems are then smoothly obtained. 
These results have potential applications in uniform controllability, the approximation of the control 
problem and state reconstruction etc. It is worthy to be investigated at length in the further research. 
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